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Abstract. In this article, we present a bifurcation and stability analysis on 
the double-diffusive convection. The main objective is to study 1) the mech- 
anism of the saddle-node bifurcation and hysteresis for the problem, 2) the 
formation, stability and transitions of the typical convection structures, and 
3) the stability of solutions. It is proved in particular that there are two differ- 
ent types of transitions: continuous and jump, which are determined explicitly 
using some physical relevant nondimensional parameters. It is also proved that 
the jump transition always leads to the existence of a saddle-node bifurcation 
and hysteresis phenomena. 



1. Introduction 

Convective motions occur in a fluid when there are density variations present. 
Double-diffusive convection is the name given to such convective motions when the 
density variations are caused by two different components which have different rates 
of diffusion. Double-diffusion was first originally discovered in the 1857 by Jevons 
[4], forgotten, and then rediscovered as an "oceanographic curiosity" a century later; 
see among others Stommel, Arons and Blanchard T2] , Veronis [IT , and Baines and 
Gill jl]. In addition to its effects on oceanic circulation, double-diffusion convection 
has wide applications to such diverse fields as growing crystals, the dynamics of 
magma chambers and convection in the sun. 

The best known double-diffusive instabilities are "salt-fingers" as discussed in 
the pioneering work by Stern [llj . These arise when hot salty water lies over cold 
fresh water of a higher density and consist of long fingers of rising and sinking water. 
A blob of hot salty water which finds itself surrounded by cold fresh water rapidly 
loses its heat while retaining its salt due to the very different rates of diffusion 
of heat and salt. The blob becomes cold and salty and hence denser than the 
surrounding fluid. This tends to make the blob sink further, drawing down more 
hot salty water from above giving rise to sinking fingers of fluid. 

The main objective of this article is to develop a bifurcation and stability theory 
for the double-diffusive convection, including 

1) existence of bifurcations/transitions, 

2) asymptotic stability of bifurcated solutions, and 

3) the structure/patterns and their stability/transitions in the physical space. 
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The analysis is based on a bifurcation theory for nonhnear partial differential 
equations and a geometric theory of two-dimensional (2D) incompressible flows, 
both developed recently by two of the authors; see respectively [71 [5] and [3] and 
the references therein. 

This bifurcation theory is centered at a new notion of bifurcation, called attractor 
bifurcation for dynamical systems, both finite dimensional and infinite dimensional. 
The main ingredients of the theory include a) the attractor bifurcation theory, b) 
steady state bifurcation for a class of nonlinear problems with even order non- 
degenerate nonlinearities, regardless of the multiplicity of the eigenvalues, and c) 
new strategies for the Lyapunov-Schmidt reduction and the center manifold reduc- 
tion procedures. The bifurcation theory has been applied to various problems from 
science and engineering, including, in particular, the Kuramoto-Sivashinshy equa- 
tion, the Cahn-Hillard equation, the Ginzburg-Landau equation, Reaction-Diffusion 
equations in Biology and Chemistry, the Benard convection and the Taylor problem 
in fluid dynamics. 

The geometric theory of 2D incompressible flows was initiated by the authors 
to study the structure and its stability and transitions of 2-D incompressible fluid 
flows in the physical spaces. This program of study consists of research in direc- 
tions: 1) the study of the structure and its transitions/evolutions of divergence- free 
vector flelds, and 2) the study of the structure and its transitions of velocity fields 
for 2-D incompressible fluid flows governed by the Navier-Stokes equations or the 
Euler equations. The study in Area 1) is more kinematic in nature, and the results 
and methods developed can naturally be applied to other problems of mathematical 
physics involving divergence-free vector fields. In fluid dynamics context, the study 
in Area 2) involves specific connections between the solutions of the Navier-Stokes 
or the Euler equations and fiow structure in the physical space. In other words, 
this area of research links the kinematics to the dynamics of fluid flows. This is 
unquestionably an important and difficult problem. Progresses have been made in 
several directions. First, a new rigorous characterization of boundary layer sepa- 
rations for 2D viscous incompressible flows is developed recently by the authors, 
in collaboration in part with Michael Ghil; see [9J and the references therein. An- 
other example in this area is the structure (e.g. rolls) in the physical space in the 
Rayleigh-Benard convection, using the structural stability theorem developed in 
Area 1) together with the application of the aforementioned bifurcation theory; see 

mm- 

In this article, we consider two-dimensional double-diffusive convections mod- 
elled by the Boussinesq equations with two diffusion equations of the temperature 
and salinity functions. In comparison to the Rayleigh-Benard convection case, and 
the steady linearized problem around the basic state for the double-diffusive con- 
vection is nonsymmetric. This leads to a much harder eigenvalue problem, and 
consequently much more involved bifurcation and stability analysis. Hence, the 
bifurcation and the flow structure are much richer. 

The central gravity of the analysis is the reduction of the problem to the center 
manifold in the first unstable eigendirections, based on an approximation formula 
for the center manifold function. The key idea is to find the approximation of the 
reduction to certain order, leading to a "nondegenerate" system with higher order 
perturbations. The full bifurcation and stability analysis are then carried out using 
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a combination of the attractor bifurcation theory and the geometric theory of 2D 
incompressible flows. 

We now address briefly the main characteristics of the bifurcation and stability 
analysis of the two-dimensional double-diffusive model presented in this article. 

First, the double-diffusive system involves four important nondimensional pa- 
rameters: the thermal Rayleigh number A, the solute Rayleigh number i], the 
Prandtl number a and the Lewis number r, defined by (|2.5p . We examine in 
this article different transition/instability regimes defined by these parameters. We 
aim to get a better understanding of the different physical mechanisms involved in 
the onset of convection. It is hoped that this will enable progress to be made in 
the theoretical understanding of the onset of double diffusive instabilities. 

From the physical point of view, it is natural to consider only the case where the 
Prandtl number a > 1. The Lewis number r measures the ratio of two diffusivities. 
From the oceanic circulation point of view, the heat diffuses about 100 times more 
rapidly than salt [TT]; hence r < 1. In this case, different regimes of stabilities 
and instabilities/transitions of the basic state can be described by regions in the 
X-rj plane (the thermal and salt Rayleigh numbers) as shown in Figure [37T] In this 
article, we focus on the regimes where 

07 

(1.1) ^<^^ = _^V(l + a-i)(l-r)-i. 

In the case where rj > rjc, transitions to periodic or aperiodic solutions are expected, 
and will be addressed elsewhere. 

Second, we show that there are two different transition regimes: continuous 
and jump, dictated by a nondimensional parameter 

(1.2) ^^^^|,V(i_^2)-l^ 

For the regime with rj < 77^ , the transition is continuous when the thermal 
Rayleigh number A crosses a critical value 

(1.3) XM = - + x'^'- 

T 4 

The rigorous result in this case is stated in Theorem 13.31 

For the regime with ry^ < rj < rjc, the transition is jump near A = Xdrj). It is 
shown also that for this case, there is a saddle- node bifurcation at A*(ry) < Ac(r?), 
together with the hysteresis feature of the transition in A* (77) < A < Ac(77). The 
rigorous result in this case is stated in Theorem 13. 4[ and schematically illustrated 
by Figure [3?3l 

Third, as an attractor, the bifurcated attractor has asymptotic stability in the 
sense that it attracts all solutions with initial data in the phase space outside 
of the stable manifold of the basic state. As Kirchgassner indicated in an 
ideal stability theorem would include all physically meaningful perturbations and 
establish the local stability of a selected class of stationary solutions, and today 
we are still far from this goal. On the other hand, fluid flows are normally time 
dependent. Therefore bifurcation analysis for steady state problems provides in 
general only partial answers to the problem, and is not enough for solving the 
stability problem. Hence it appears that the right notion of asymptotic stability 
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should be best described by the attractor near, but excluding, the trivial state. It 
is one of our main motivations for introducing the attractor bifurcation. 

Fourth, another important aspect of the study is to classify the structure/pattern 
of the solutions after the bifurcation. A natural tool to attack this problem is the 
structural stability of the solutions in the physical space. Thanks to the aforemen- 
tioned geometric theory of 2D incompressible flows, the structure and its transitions 
of the convection states in the physical space is analyzed, leading in particular to a 
rigorous justification of the roll structures. More patterns and structures associated 
with the double diffusive models will be studied elsewhere. 

Fifth, for mathematical completeness and for applications to other physical 
problems, the bifurcation analysis is also carried out for the case where the Lewis 
number r > 1. In this case, only the continuous transition is present. The rigorous 
result in this case is stated in Theorem 13.51 

This article is organized as follows. The basic governing equations are given in 
Section 2, and the main theorems are stated in Section 3. The remaining sections 
are devoted to the proof of the main theorems, with Section 4 on a recapitulation of 
the attractor bifurcation theory and the geometric theory of incompressible flows. 
Section 5 on eigenvalue problems. Section 6 on center manifold reductions, and 
Section 7 on the completion of the proofs. 



2. Equations and Set-up 



2.1. Boussinesq equations. In this paper, we consider the double-diffusive con- 
vection problem in a two-dimensional (2D) domain x (0, h) C (ft, > 0) with 
coordinates denoted by (x, z). The Boussinesq equations, which govern the motion 
and states of the fluid flow, are as follows; see Veronis ^14j : 

^ + {U- V)U = - — (Vp + pqe) + z/Ai7, 
ot po 



(2.1) 



dT 

— + ({/• v)r - ktAt, 

ot 

— + {U-V)S^KsAS, 
ot 

div U = 0, 



where U = {u, w) is the velocity function, T is the temperature function, S is the 
solute concentration, P is the pressure, g is the gravity constant, e = (0, 1) is the 
unit vector in the z-direction, the constant i/ > is the kinematic viscosity, the 
constant kt > is the thermal diffusity, the constant K5 > is the solute diffusity, 
the constant po > is the fluid density at the lower surface z — 0, and p is the fluid 
density given by the following equation of state 



(2.2) 



p = Po[l-a{T-To) + b{S-So)]. 



Here a and b are assumed to be positive constants. Moreover, the lower boundary 
{z — 0) is maintained at a constant temperature Tq and a constant solute concentra- 
tion 5*0, while the upper boundary (z = h) is maintained at a constant temperature 
Ti and a constant solute concentration 5i , where Tq > Ti and > 5*1 . The trivial 
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steady state solution of (|2.ip - ()2.2|) is given by 



(2.3) 



^ [/" = 0, 



where po is a constant. To make the equations non-dimensional, we consider the 
perturbation of the solution from the trivial convection state: 



U" = U-U°, y„^y_j,o^ 

Then we set 

{x, z) ~ h{x' , z'), t = h^t'/KT, 

U" = ktU'/K T" = (To - Ti)T', 

S" = [So - Si)S', p" - poi'KTPV/i'. 

Omitting the primes, the equations (|2.ip can be written as 

^ = a(AC/ - Vp) + (j{\T - T]S)e - {U ■ V)f/, 
dT 

— = AT + w-{U-V)T, 
tAS + w-{U- W)S, 



(2.4) 



as 

'dt 
divU = 0, 



for {x,z) in the non-dimensional domain f2 = x (0, 1), where U — {u,w), and 
the positive nondimensional parameters used above are given by 



(2.5) 



A = 



ag{To ~ Ti)h^ 
bg{So ~ Si)h^ 



kt 



the thermal Rayleigh number, 

the salinity Rayleigh number, 
the Prandtl number, 
the Lewis number. 



We consider periodic boundary condition in the x-direction 

(2.6) {U,T,S){x,z,t) {U,T,S){x + 2kTT/a,z,t) V/s G Z. 

At the top and bottom boundaries, we impose the free- free boundary conditions 
given by 



(2.7) 



3u 

(T,S,w)^0, T^ = at z = 0,l. 
oz 
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It's natural to put the constraint 
(2.8) / udxdz = 0, 



for the problem (|2.4p - (|2.7p . It is easy to see from the following computation that 
(|2.4p is invariant under this constraint: 

Q pi pin j a pi pin j a 

— / udxdz = / / (j{uxx + Uzz -Px) - (uux + wuz)dxdz 
OtJo Jo Jo Jo 

1 plrr/a 

o-{ux -p) 11^=0''^" dz+ a{uz \lZ:l)dx 

Jo 

1 pItt/o 

/ {ux + Wz)udxdz 
Jo 
=0. 

The initial value conditions are given by 
(2.9) {U,T,S) ^ {U,f,S) at i = 0. 

2.2. Functional setting. Let 

H ={{U,T,S) £ L^{n)^ I dWU ^0,w\z^o,i^O, I udxdz = 0, 

Jn 

u is 27r/a-periodic in a;-direction }, 
V ={([/, T, S) e H^iVlf r\H I ([/, r, S) is 27r/a-periodic in x-direction, 
T \z=o,i= S 1^=0,1= 0}, 
Hi ^VDH^in). 
Let G . Hi ^ H and L^jj = —A — B\,^ : Hi ^ H he defined by 

G(V') - i-pm ■ v)u],~{u ■ v)r, ~{u ■ v)5), 

Bxrj-ip = i-P[a{XT - r]S)e], -w, -w), 

for any ip = {U,T, S) E Hi. Here P is the Leray projection to fields. Then the 
Boussinesq equations (|2.4p - (l2.8p can be written in the following operator form 



(2.10) ^ = La^^ + G(^), ^ = iU,T,S). 

3. Main Results 

3.1. Definition of attractor bifurcation. In order to state the main theorems of 
this article, we proceed with the definition of attractor bifurcation, first introduced 
by two of the authors in [SI . 

Let H and Hi be two Hilbert spaces, and Hi ^ H he a, dense and compact 
inclusion. We consider the following nonlinear evolution equations 

(3.1) ^^ = L,u + Giu,X), 

u{0) = Mo, 
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where u : [0,oo) H is the unknown function, A G M is the system parameter, 
and Lx : Hi H are parameterized hnear completely continuous fields depending 
continuously on A G R-'^, which satisfy 

{— L\ = A + B\ a sectorial operator, 
A : Hi ^ H a linear homeomorphism, 
Bx : Hi H parameterized linear compact operators. 

It is easy to see [3] that Lx generates an analytic semi-group {e*^^}t>o. Then 
we can define fractional power operators {~Lx)^ for any < M < 1 with domain 
-ff^ = D{{-LxY) such that H^^ C H^^ if > /i2, and Hq = H. 

Furthermore, we assume that the nonlinear terms G'(-,A) : iJ^ — > H for some 
1 > /i > are a family of parameterized C" bounded operators {r >1) continuously 
depending on the parameter A € M^, such that 

(3.3) G{u,\) = o{\\u\\h^X VAeRi. 

In this paper, we are interested in the sectorial operator —Lx = A + Bx such 
that there exist an eigenvalue sequence {pk\ C C"'^ and an eigenvector sequence 
{efe,/ife} C Hi oiA: 

'Azk=PkZk, Zk=ek+ihk, 

(3.4) I Repk ^ oo {k ^ oo), 

^ |Impfe/(a + Rcpfe)| < c, 

for some a, c > 0, such that {e^, hk} is a basis of H. 

Condition (|3.4p implies that A is a sectorial operator. For the operator Bx : 
Hi ^ H, we also assume that there is a constant < 6* < 1 such that 

(3.5) Bx: Hg — > H bounded, V A e 



Under conditions (|3.4p and (|3.5p . the operator —Lx — A+Bx is a sectorial operator. 

Let {Sx{t)}t>o be an operator semi-group generated by the equation (|3.ip . Then 
the solution of (|3.ip can be expressed as 

^it,^Q) = Sxit)ibo, t>0. 

Definition 3.1. A set T, C H is called an invariant set of iS.l]) if S{t)H = E 
for any t > 0. An invariant set T, <Z H of USA]) is said to he an attractor if S is 
compact, and there exists a neighborhood W <Z H of Y, such that for any t/jq £ W 
we have 

lim disti/(-!/'(i, V'o), S) — 0. 

t — >oo 

Definition 3.2. (1) We say that the solution to equation bifurcates from 
(-0, A) — (0, Ao) to an invariant set fix, if there exists a sequence of invariant 
sets {f^A„} of liS.l]} such that ^ f^A„, md 

lim A„ = Ao, 

n — *oo 

lim max = 0. 

n— ►CO x€fl\„ 

(2) // the invariant sets fix are attractors of ^3.1]) . then the bifurcation is called 
attractor bifurcation. 
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3.2. Main theorems. We now consider the double diffusive equations ()2.4p . In 
this article, we always consider the case where the parameters A and 77 satisfying 

27 r-^){l-r)-\ 



(3.6) 



A« A, = 



T] 27 



T 4 

First we consider a more physically relevant diffusive regime where the thermal 
Prandtl number a is bigger than 1, and the Lewis number t is less than 1: 

(3.7) o- > 1 > r, ^ Tr^/2. 

Here the condition on a defines the aspect ratio of the domain. In this case, we 
consider two straight lines in the A — r/ parameter plane as shown in Figure 13.11 




Tic. Tic 



Figure 3.1. 



A = A,(7y), 



W ^ ^ ^ 27 4 



(3.8) 
where 

(3.9) 



Also shown in Figure 13.11 are two critical values for rj given by 

= |7rV2(l + a-i)(l - r)-\ = |^V^(1 - r". 

The following two main theorems study the transitions/bifurcation of the double- 
diffusive model near the line Li for ry < rye- 
Theorem 3.3. Assume that the condition ^3. 7[ ) holds true, andrj < 77^ ~ ^7r^T'^(l- 
T^)^^. r/ien f/ie following assertions for the problem ^2.4^ - (K^) hold true. 

(1) If X < Ac; the steady state {U,T,S) = is locally asymptotically stable for 
the problem. 

(2) The solutions bifurcate from (([/, T, S), A) = (0, Ac) to an attractor Sa for 
X > Ac, which is homeomorphic to , and consists of steady state solutions 
of the problem. 
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(3) For any "tpo — (U, T, S) G H\T, there exists a time to > such that for any 
t > to, the vector field U{t, ^o) is topologically equivalent to the structure as 
shown in Figure [X31 . where ip = (C/(t, i/jo), T(i, i/iq); 'S'(t, V'o)) is a solution 
of ^2.4^ - (273\) with ^3.7^ - (37B\) . T is the stable manifold of the trivial solution 
{U,T, S) — Q with co-dimension 2 in H. 



H 




Figure 3.2. If r > 1 or 77 < 77^, the equations bifurcate from 
(0, Ac) to an attractor for A > Ac. 



Theorem 3.4. Assume that the condition \3. ?| ) holds true, and rj^ > rj > rj^^ = 
^7r^r'^(l — T^)~^. Then there exists a saddle-node bifurcation point Aq (Xo < ^c) 
for the equations, such that the following statements for the problems ^K^-^KW 
hold true. 



(1) At X = Xo, there is an invariant set Sq = Sa,, with ^ Sq- 

(2) For A < Aq, there is no invariant set near Eq. 

(3) For Aq < A < Ac, there are two branches of invariant sets Yi\ and S^, and 
S| extends to X > Xc and near Ac as well. 

(4) For each X > Aq, is an attractor with c?ist(E^,0) > 0. 

(5) For Aq < A < Ac, 

(a) Tt\ is a repeller with ^ Tt\, and 

(b) when X is near Ac, Tt\ is homeomorphic to , consisting of steady 
states. 
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Figure 3.3. If r < 1 and f]c^ < t] < t]c, the equations have a 
saddle-node bifurcation for A < Ac- 



We now consider the diffusive parameter regime where a > 1, t > 1, = /2 
and (7 7^ r. In this case, two lines are shown in Figure 13.41 The foUowin theorem 
provides bifurcation when A crosses the line Li. 

Theorem 3.5. Assume that a > 1, t > 1, — 7r^/2, a ^ t and i3. 6]) hold, then 
for any r/ > 0, the following assertions for the problem ^2.4^ - [K^) hold true. 

(1) If X < Ac, the steady state ([/, T, S) — is locally asymptotically stable for 
the problem. 

(2) The solutions bifurcate from {{U,T, S), X) = (0, Ac) to an attractor for 
X > Xc, which is homologically equivalent to , and consists of steady state 
solutions of the problem. 

(3) For any -00 — {U,T,S) S H\T, there exists a time tQ>0 such that for any 
t > to, the vector field U{t,ipo) is topologically equivalent to the structure as 
shown in Figure fXSI . where — {U(t,ipQ),T(t,ipQ), Slt^ipo)) is a solution 
of i2.4^ - fK^) . r is the stable manifold of the trivial solution {U,T, S) — Q 
with co-dimension 2 in H. 




Figure 3.4. 
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Figure 3.5. 
4. Preliminaries 



4.1. Attractor bifurcation theory. Consider p.ip satisfying (|3.2p and (|3.3p . 

We start with the Principle of Exchange of Stabihties (PES). Let the eigenvalues 
(counting the multiplicity) of Lx be given by 

/3i(A),/32(A),-- - ,/3fe(A),---eC. 

Suppose that 

' < if A < Ao 

(4.1) Re(3i{\) <! = if A = Ao (1 < i < m) 

> if A > Ao 



(4.2) RePj{Xo) < 0. 

Let the eigenspace of L\ at Ao be 



V m + 1 < J- 



En 



[j \J{u,veHi\ {Lxo - (3j{Xo))^w ^ 0,w ^ u + iv}. 

l<j<rn k=l 

It is known that dimi?o = m. 

Theorem 4.1 (T. Ma and S. Wang [SllZj). Assume that the conditions ^^-(K^ 
and j/^. j[ )- |7^ hold true, and u = is locally asymptotically stable for i3. 1\) at 
A = Ao. Then the following assertions hold true. 

(1) \3.1\) bifurcates from (it, A) = (0, Ao) to attractors Y^x, having the same 
homology as S"^~^ , for X > Xq, with m — 1 < dimY^x < m, which is 
connected as m > 1; 

(2) For any ux G Yx, ux can be expressed as 

ux = vx+o{\\vx\\hi), I'xeEo; 

(3) There is an open set U d H with (z U such that the attractor Yx bifurcated 
from (0, Ao) attracts U\T in H , where F is the stable manifold of u = Q with 
co-dimension m. 

In the case where to = 2, the bifurcated attractor can be further classified. 
Consider a two-dimensional system as follows: 

dx 



(4.3) 



dt 



— (3(X)x — g{x,X), xE 



Here /3(A) is a continuous function of A satisfying 







^<0 


if A < Ao 


(4.4) 


/3(A) < 


= 


if A = Ao 






l>0 


if A > Ao 
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and 

r g{x,X)=gk{x,X)+o{\x\''), 

(4.5) < gk{-, A) is a k-multilinear field, 

[ Cilxl'^+i << gk{x,X),x>< C2\x\''+\ 

for some integer k ^ 2m + 1 > 3, and some constants C2 > Ci > 0. 

The following theorem was proved in JT] , which shows that under conditions 
and (|4.5p . the system (|4.3p bifurcates to an S'^-attractor. 

Theorem 4.2. Let the conditions {4-4^ and \4.5^ hold true. Then the solution to 
the system {4 .3^ bifurcates from (x, A) — (0, Aq) to an attractor Ea for A > Aq, 
which is homeomorphic to . Moreover, one and only one of the following is true. 

(1) Ea is a periodic orbit, 

(2) Tix consists of only singular points, or 

(3) Sa contains at most 2{k + 1) = 4(m + 1) singular points, and has AN + n 
(N + n > 1) singular points, 2N of which are saddle points, 2N of which 
are stable node points (possibly degenerate), andn of which have index zero. 

4.2. Center manifold reduction. A crucial ingredient for the proof of the main 
theorems using the above attractor bifurcation theorems is an approximation for- 
mula for center manifold functions derived in [7] . 
Let Hi and H be decomposed into 



(4.6) 

for A near Aq G M^, where E^, E2 are invariant subspaces of L\, such that 



H = E^®E^, 



dim_E;^ < 00, 
E^ - E^, 



E2 = closure of E2 in H. 



In addition, La can be decomposed into L\ = C'^ ® C2 such that for any A near Aq, 
(4.7) 




where all eigenvalues of £3 possess negative real parts, and the eigenvalues of C\ 
possess nonnegative real parts at A = Aq. 

Thus, for A near Aq, equation p.ip can be written as 

5 -/:^a; + Gi(a;,y,A), 



(4.8) 



dt 

^=£^y + G2{x,y,X), 



where u = x + y & Hi, x & E^, y e E^, Gi{x, y, A) = PiG{u, A), and : H ^ E^ 
are canonical projections. Furthermore, let 

i?2 (/i) = closure of E2 in H^ , 

with /i < 1 given by (|3.3p . 
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By the classical center manifold theorem (see among others t3l[I3)i there exists 
a neighborhood of Ao given by |A — Ao| < S for some 6 > 0, a neighborhood B\ C 
of a; = 0, and a center manifold function <!>(•, X) : B\ E2{9), called the center 
manifold function, depending continuously on A. Then to investigate the dynamic 
bifurcation of (|3.ip it suffices to consider the finite dimensional system as follows 

(4.9) ^ = £^x + gi(x,a>;,(x),A), xeB^cE^. 

Let the nonlinear operator G be in the following form 

(4.10) G{u,X)^Gk{u,\)+o{\\u\\''), 
for some integer k > 2. Here Gk is a fc-multilinear operator 

Gk:HiX---xHi^H, 
Gfe(u, A) = Gk{u,- ■ ■ ,u, A). 

Theorem 4.3. [Tj Under the conditions ^-6^ , J^. 7p and jO[ ), the center manifold 
function $(x, A) can be expressed as 

(4.11) ^x, L) = {-C^)-^P2Gk{x, A) + o{\\xt) + 0(|Re/3| ||xf ), 

where £3 *^ |.^.7| j, P2 '■ H E2 the canonical projection, x G E^ , and 

(3 = (/3i(A), • • • ,/3m(A)) t/ie eigenvectors o/£^. 

Remark 4.4. Suppose that {ej}j, the (generalized) eigenvectors of L\, form a 
basis of H with the dual basis {e*}j such that 

= if i^j, 
^0 if i-j. 



(ei,e*)ij 



Then, we have 



m 



i=l 



Hence, near A = Aq, p2Gk{x, A) can be expressed as follows. 

OG 

(4.12) P2G{x,X)^ E Gi(x,A)e,+o(||x||'=), 

j=m+l 

where 

(a; , A) = E a^-^ . . .^-^ a^j, • • • Xj^ , 
= (Gkiej,,- ■ ■ ,ej,, A),e*)///(ej,e*). 



eigenvalues /?i(A), • • • , (3„i{X) satisfy 



In many applications, the coefficients can be computed, and the first m 



|Re/3(Ao) 



E(Re/3,(Ao))2 = 0. 
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Hence (|4.11|) and (I4.12p give an explicit formula for the first approximation of the 
center manifold functions. 

4.3. Structural stability theorems. In this subsection, we recall some results 
on structural stability for 2D divergence- free vector fields developed in which 
are crucial to study the asymptotic structure in the physical space of the bifurcated 
solutions of the double-diffusive problem. 

Let C"'(17,M2) be the space of all C"^ (r > 1) vector fields on f7 = x (0, 1), 
which are periodic in x direction with period 2Tr/a, let D^{il, R^) be the space of all 
C divergence- free vector fields on O = R^ x (0, 1), which are periodic in x direction 
with period 27r/a, and with no normal flow condition in z-direction: 

D'^ifl^R^) = {v= iu,w){x,z) e C"'(r2,R2) I w = at z = 0, 1 } . 

Furthermore, we let 

Bl^{n,R^) = {w e D''{n,R^) I w = at z ==0,1 } , 

f u = at z = 1 

w — — = at z = 1 
oz ■> 

Definition 4.5. Two vector fields Vi,V2 G C""(f2,R^) are called topologically equiv- 
alent if there exists a homeomorphism ofip:Q-^fl, which takes the orbits of vi to 
orbits of V2 and preserves their orientation. 

Definition 4.6. Let X = D''{n,M?) or X = BliVL^m?). A vector field vq G X is 
called structurally stable in X if there exists a neighborhood U G X of Vq such that 
for any v Cz U , Vq and v are topologically equivalent. 

Let V E D^{n, R^). We recall next some basic facts and definitions on divergence- 
free vector fields. 

(1) A point p G 57 is called a singular point of v if v{p) = 0; a singular point p 
of V is called non-degenerate if the Jacobian matrix Dv{p) is invertible; v 
is called regular if all singular points of v are non-degenerate. 

(2) An interior non-degenerate singular point of v can be either a center or a 
saddle, and a non-degenerate boundary singularity must be a saddle. 

(3) Saddles of v must be connected to saddles. An interior saddle p G is 
called self-connected if p is connected only to itself, i.e., p occurs in a graph 
whose topological form is that of the number 8. 

Theorem 4.7. For vector fields satisfying free-free boundary conditions, we set 

du 



Q at z^Q.l 
oz 



udxdz = 



B5(17,R2) = <^ i> g B5(f7,R2) 

I " Jo. 

Then v G Bl^iVL.M?) (resp.v G Bl{n,M?)) is structurally stable in B5(f7,R2) (resp. 

in i32(f7,R^)y' if and only if 

1) V is regular; 

2) all interior saddle points of v are self-connected; and 

3) each boundary saddle of v is connected to boundary saddles on the same 
connected component of dQ ( resp. each boundary saddle of v is connected 
to boundary saddles not necessarily on the same connected component) . 
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5. Eigenvalue problem 



(5.1) 



In order to apply the center manifold theory to reduce the bifurcation problems, 
we shall analyze the following eigenvalue problem for the linearized equations of 

(El-dlii). 

( a{AU - Vp) + {aXT - aTjS)e = PU, 
AT + w^ (3T, 
tAS + w^pS, 
divU = 0, 
Ou 

-Q^ \z=o,i— w |z=oa= T l^^g 1= S \z=o,i— 0. 

We prove couples of lemmas to show that the operators —L\f^ are sectorial operators 
when the parameters are properly chosen. In order to get the precise form of the 
center manifold reduction, the eigenspaces are analyzed in detail in this section. 

5.1. Eigenvalues. We shall use the method of separation of variables to deal with 
problem (|5.1[) . Since tp ~ iU,T, S) is periodic in x-direction with period 2TT/a, we 
expand the fields in Fourier series as 

oo 

(5.2) i'{x,z)= M^y^^- 

j=-oo 

Plugging (|5.2p into (|5.ip . we obtain the following system of ordinary differential 
equations 

' DjUj — ijapj = (J~^ I3uj, 



(5.3) 



- Wi 



\Tj - riSj 



(J ^PWj, 



rDjSj + Wj = pS- 



J' 



I J auj 



- Wa 



0, 



^ u'j l^^o i— Wj l^^oa— Tj 1^^0,1— Sj \z=o,i— 



for j G Z, where ' = d/dz, Dj — d /dz —j a . We reduce 
for Wj{z): 

(5.4) 



to a single equation 



{{rD,-p)(D,-(3){D,-a-^(3)D, 

+ pa^[\{TD, - /3) - r^{D, ~ P)]}w, = 0, 



(5.5) 



Wi 



w 



(4) 



W 



(6) 



at z = 0,l, 



for j G Z. Thanks to (|5.5p . Wj can be expanded in a Fourier sine series 

oo 

(5.6) Wj (z) = Wji sin lirz 

1=1 

for j G Z. Substituting (|5.6p into (|5.4p . we see that the eigenvalues f3 of the problem 
(|5.ip satisfy the following cubic equations 

(5.7) (3' + (<7 + T + 1)7^/32 + [{a + T + nT)^% - afa^^J^^X - ry)]/? + ar-fl 

+ (jfa'^irj - tA) = 0, 
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for j e Z and I £ N, where — j^a^ + P-k"^. 

For the sake of convenience to analyze the distribution of the eigenvalues, we 
make the following definitions. 

Definition 5.1. For fixed parameters a, t, rj and \, let 

(1) = /33 + (a + r + + (a + T + aT)7j,/3 + aT^%, 

(2) h,i{P) = [afa^j7^\X ~ 77)]/3 - afa^T^ ~ rX), 

(3) fjiW)^9dP)~hjiiP), 

(4) ^, = f7r4r2(l + c7-i)(l-r)-i, 

(5) ^e, =f7rV3(l-r2)-i, 

(6) A, = f + f ^^ 

(7) (3jii,(3ji2 and P J 13 be the zeros of fji with Re{(3jii) > Re{(3ji2) > Re{(3ju). 
In the following discussions, we shall focus on the following diffusive regime: 

(5.8) cr>I>r>0, 7r^/2, 77 < 77c and A w Ac. 

Lemma 5.2. (1) Under the assumption \5.l^] . has three simple real ze- 

ros. 

(2) If T > 1, Pill — is a simple zero o//ii(/3) for A — Ac. 
Proof. Since A w Ac, it suffices to prove this statement for A = Ac. In this case, 
/ii(/3) + [STT^a + T + I)/2]/32 + [97r\a + t + ar)/A - a(Ac - r;)/3]/3 
=/(/3)/3. 

Since rj < rjc or t > f, the constant term of is nonzero. Hence Pm = is 
a simple zero of fn. Moreover, for condition (|5.8|) . the quadratic discriminant of 
/(/?) is 97r4(cr + I - t)2/4 + 4(Tr/(I - r)/(3r) > . This implies fn has three simple 
real zeros. □ 

We summarize the following important lemma about the distribution of the zeros 

of fjl- 

Lemma 5.3. Assume that either 

1) rj < rjc with condition fi5.8\) or 

2) 77 > with T >l, 

then 



< 


if \< Ac, 


= 


i/ A = Ac, 


> 


if \> Ac, 



(5.9) PiiiiX) 

(5.10) RePMX) < for (j, I, k) ^ (I, I, I). 
Proof. Let 

f{p) = + a2/3' + aip + ao 
be a monic real coefficient polynomial of degree 3. The following cases are apparent. 
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(1) /3 = is a zero of /(/3) if and only if ao — 0. 

(2) /3 — bi,~bi,a {a,b E M) are zeros of f{/3) if and only if ai = b^, 02 = —a 
and ao = — a6^. 

Case (2) is equivalent to 

ai > and aia2 — ag. 

By the above observation, we prove the lemma in several steps as follows. 

Step 1. It's easy to see that (3 = —70;, —<^loi, —tJoi zeros of /o;. When 

j 7^ 0, /? = is a zero of fji if and only if 

<^TJji + crfa'^iv - tX) = 0, 

which is equivalent to 

(5.11) A=^ + ^. 



For fixed rj, minimizing the right hand side of (|5.1ip . we obtain that Ac = ^ + ^tt* 
is the global minimum of A in Case (1), provided — ^ and (j, I) ~ (1, 1). 
Step 2. For Case (2), we obtain in the same fashion as above that 



(u + 1) ' apa^ ' 



(5.12) A = + ^(a + r)(r + 1). 



For a fixed 77, minimizing the right hand side of (|5.12p . we obtain that 
A.=|^.+ f.^(l + .-V)(l + .) 

is the global minimum of A in the case (2), provided ^ and (j, I) = (1, 1). 

Step 3. As introduced before, A — Ac(r/) and A = Xc^{rj) define two straight 

lines in the X — rj plane, shown in Figures l3. II and 13.41 

If T < 1, the intersection of the two lines is 

27 27 
(Ac., 77c) = (^ttMI + ra-i)(l - r), -j^'t\1 + a-'){l - r)-'). 

As shown in Figure [3711 Xc{r]) < Aci (??) for rj < rjc- If r > 1, Xc{r]) < Aci(?7) for 
7/ > 0; see Figure It is easy then to see that in either case, (15. 9p and (|5.10p 
hold true. The proof is complete. □ 

Remark 5.4. (1) This lemma works for the 3D double-diffusive problem as 
well. 

(2) In case 1), if 77 > 77c then Pm = /3ii2 are complex numbers for A ~ A^. 

(3) The distribution of the zeros of fji was first analyzed by Veronis 14 . The 
results are scattered in different papers. To make this paper more self- 
contained, the authors think it's good to summarize it here and give a clear 
proof. 

To check that the operators —L\^ satisfy condition p.4p . we prove the following 
lemma. 

Lemma 5.5. (1) Only finitely many numbers of the zeros of fji{P) have nonzero 
imaginary parts for (j, Z) G Z x N. 
(2) Pjik -00 as f + P 00. 
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Proof. Since fji — gji — hji , /3 is a zero of fji (/?) if and only if /3 satisfies the equation 

(5.13) g,i{(3) = h,i{p). 
Plugging /3 = jjiP* into (|5.13p . we obtain 

(5.14) (/3* + l)(/3* + r)(/3* + a) = z?,,[(A - 77)/?* - {rj - rA)], 

where i?,/ — pa^a/^^. Since lim dji ~ 0, the roots of (|5.14[) must be negative 

real numbers near the interval [—a, — r] when (j^ + P) is large. This completes the 
proof. □ 

5.2. Eigenvectors. Let's make some observations to analyze the spectrum. Since 
9df3) = (/3 + 7|/)(/3 + rj]Mf3 + 0^%) and h,i = crj^a^iJi^X - v)f3 - (v - rX)j],], 
it's easy to check that (3 = —'Jji or /3 = —tj^i is a zero of fjiiP) if and only if j — 0. 
In the case of j = 0, the zeros of fji are —Jji, ^nd —<tjji. The corresponding 

eigenvectors are 

(5.15) Voi(2^>2) = (0,0,sin?7rz,0)*, 

ilj'^i{x,z) = (0, 0, 0, sinZTTz)*, 

lpQl{x,z) = (cos ^TTZ, 0,0,0)*. 

To analyze the structures of the eigenspaces of problem (|5.ip . we make the following 
definitions. 

Definition 5.6. For j ^ 0, we define 
In 

(l)]i{x,z) = ( — cos jax cos Zttz, sin jax sin /ttz, 0, 0)*, 
ja 

(t)^i{x,z) = (0, 0, sin jax sin Zttz, 0)*, (j)'ji{x,z) — (0, 0, 0, sin jax sin Zttz)*, 
Itt 
ja 

(t)^i{x,z) = (0, 0, COS jax sin Zttz, 0)*, (j)^i{x,z) = (0, 0, 0, cos jax sin /ttz)*. 
for each I G N. 

Lemma 5.7. // j 7^ and [3 is a zero of fji, then we have the followings. 
(1) The eigenvector corresponding to (3 in the complexified space is 

(5.16) ^^,(x,z) = e''^"'"(— cos Itt z, sin Itt z,Ai{P) sin Itt z,A2(p) sin In zY, 

^ ja 

1 1 

where ^i(/3) — — j- and ^2(/5) 



(f)'^i{x,z) — ( — ^ sin jax cos ^TTz, cos jax sin ^7r2;, 0, 0)*, 



(2) If (3 is a real number, the corresponding eigenvectors are given by 

(5.17) Vf/' = 0]/ + A,iP)(f>% + A2(/3)4 and 

V'f/' = 4+Ai(/3)4+A2(/3)4. 
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(3) // Im{(3) ^ 0, the generalized eigenvectors corresponding to {3 and (3 are 
5.18 

= 4>]i + RiCmi + R2Cm% + /i(/5)0f, + hCmi and 

w/iere = i?e(yli(/3)), Im{Ai{l3)), i?2(/3) = i?e(v42(/3)) a«rf 

UP) - /m(^2(/3)). 

The proof of Lemma ISTTl follows from a direct calculation, and we shall omit the 
details. 

Definition 5.8. (1) If j = and / e N, we define 

Eqi = span{ipli{x,z),tpli{x,z),ip^i{x,z)}. 

(2) For j £ N, we define Eji — span {</>]; (a;, z), (/'|;(a;, z), (fiji{x, z)}, Ej^ — span 
{(j)'^i{x, z),<j>^i{x, z), (x, z)} and E^i = Eji ® E^. 

(3) For J € U N, I S we define Ef .^ be the eigenspace spanned by the 
eigenvectors and the generalized eigenvectors corresponding to the zeros of 
fu- 
ll is easy to see that the completion of ©!^oi=i-^j' ™ H-norm is H. Hence 

the following theorem shows that the eigenvectors and the generalized eigenvectors 
corresponding to the zeros of {fji}^Q form a basis of H. 

Theorem 5.9. Under the assumption h5.^] . we have 

1) Ef^^ = Eji for j G {0} U N, Z G N; and 

2) L\^\Eji is strictly negative definite for each {j,l) G Z x N when A < Ac. 

Proof. We proceed in two steps. 

Step 1. To prove Assertion 1), it is enough to show that dim Ef^^ > dim Eji. 
The case of j = follows from (|5.15p . When j ^ 0, for a fixed I, we examine all 
the cases as follows. 

1. If Pi > (32 > are distinct zeros of fji{P), by (|5.17p . we have dim Ef.^ > 6 = 
dim Eji . 

2. If /3i /32 G C\M and are distinct zeros of fji{P), by ((STT)) and ((5l8l) . we 
have dim Ej.^ > 6 = dim Eji . 

3. The first derivative of fji is /j, = 3/3^ + 2{a + T + l)7|;/3 + [(cr + r + o-t)7^, - 
apa'^^Ji^{\ — 7^)]. The quadratic discriminant of /j, is 

^{\[{^ - rf + (a - \f + (r - \f\^% + ■io2''o?l-^^{\ n)\ > 0, 
since A w Ac > 77. It follows that fji{P) cannot have zeros of multiplicity 3. 
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4. If l3i (32 — are zeros of fji, direct computation shows 



(5.19) 



- fi2) "^'^ 



(/32 - /3i) "^'^ 



+ [ 



Note that 



afa^X{T- 1). 



Hence 



/33 = 



(aAj^a^)(r-l) 
(/3i+7|,)(/32+7|z) 



We pick up u = 'm-'P'ji to be the generaUzed eigenvector corresponding to /?3 in E^i , 



where m is some smah constant. It's easy to check that Ej^ 
and L\jj\Eji can be represented by matrix 

/a 



span{V'f/'\ V'fr\w} 



(5.20) 





P2 

\ 



' ^- (/3i-/32) I 



.-2 „ 2 



7?, (/32-/3i) 



/33 



in the basis {ipji' ^i^ji' i"^}- The same argument works for Ej^ as well. 

Step 2. It's easy to check that Ej^ is orthogonal to Eji for e N x N and 
Ej^i^ is orthogonal to £^^2^2 for {ji,li) ^ (^2,^2)- Lemma [Ql together with Stepl 
imply that Lxrj\Eji is strictly negative definite when A < Ac. This completes the 
proof. □ 

Lemma 15.31 and Lemma 15.51 together with Theorem 15.91 imply the following the- 



Theorem 5.10. Under assumption i5.8\) . —Lxri is a sectorial operator. 
Remark 5.11. (1) Since dim Eji = 3 or 6 which is finite, there exists a vector 



^'f;'' e Eji such that 



"31 



= for (/3,fc) 7^ (/3*,fc*), 
7^0 for (/3,fc) = (/3*,fc*), 



where /3 and /3* are zeros of fji and fc, k* — 1, 2. 
(2) Note that E'j; is orthogonal to E'^i for (j, /) g N x N and E'^iZ^ is orthogonal 

to Ej^i^ for (ji, ^i) 7^ (j2, ^2), hence we conclude that < *f;fi', V'l.jf >-?/= 
for {ji,li,(3,ki) ^ (j2,^2,/3*,fc2). 
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(3) If j = 0, we pick up 'i'^i ~ il^^i- For (j, I) = (1, 1), we pick up 

(5.21) = 0ii+Ci02i+C2(/.?i, and = <^ti + C^i0ii + ^s^^i, 
where 

(5.22) C,= , and C2 = --— 

/3iii + 7ii All + T7fi 

(4) Lemma 15.21 and Theorem 15.91 show that the multiphcity of the eigenvalue 
/3iii(A) is two and the corresponding eigenvectors are ■(/'fi"'^ and -i/'n"'^- 

(5) For ij] e -^3/4 C i/, by Sobolev inequahty, 

/"l 1*277 1 a 

where C is some constant. Hence, G{ip) — o{\tp\H3/^)- 
6. Center manifold reduction 



We are now in a position to reduce equations of (|2.4|) - (|2.8p to the center manifold. 
We would hke to fix 77 < 77c , and let A « Ac be the bifurcation parameter. For any 
ip = {U, r, S) e H, we have 

00 3 
j=0,i=l fc=l 

Since /^m is the first eigenvalue, V'n"'^ and V'n"'^ are the first eigenvectors. The 
reduced equations are given by 



'^''"^ /3iii(A)xni+ ,^,,,.1 <G(7/.,V),*fr''>if, 



(6.1) 



< vii ,y^ii >-ff 

Here for = (C/i, Ti, 5i), ^2 = (?72, ^2, ^2) and V>3 - {U3, T3, ^3), 

G(V'i, ^2) = -{P{Ui ■ V)U2, {Ui ■ V)T2, {Ui ■ V)S2Y and 

< G{^bi,^2), V'3 >H= - f[< {Ui ■ V)C/2, C/3 >R2 +([/i • V)T2T3 + (C/i • V)5253]da:dz, 

where P is the Leray projection to fields. 

Let the center manifold function be denoted by 

(6.2) '^^ J2 K''(^iii,yiii)V'f/V<i>f/'(a;ni,yiii)V'f/'). 

Note that for any ipi ^ Hi{i — 1, 2, 3), 

< G(V'l,l/'2),'(/'2 >H^ 0, 

< G(V'i,i/'2),i/'3 >H= - < G(7/'i, V'3), V'2 >h; 

and for k=l,2, 

< G(7Ai, <"■'=), ^fi"^'^ >H=0. 

Then by = a^iiiV'ii"'"^ + yiiiV'ii"'^ + 'J'j we have 
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(6.3) <G(^,^),vI,^i"'i >H=<G(V'f^'^ vl/^r^' >ff 2/?ii 

-<G(<"'\vI/fi"'i),$>^xin 
-<G(<"-^vI/fi"^l),$>^ 2/111 

+ <G(<i>,<-^'^),vI/?{-'^ >Hyiii 
+ <G(<i>,<i>),vI'fi"-i >H, 



(6.4) < G(V^,^),*?i"'^ >H= < G(V'fr'\ <--!), '2 >ff x?ii 

+ < G(<"-^<"'^),v^'?--2 >ff xnmu 

~ < G{i;^,r'\'^ir-^),'^ >H ^111 

+ < G($,V'fr'^),*?i"^' >HXni 
+ < G($,$),*fi"'^ >H, 

By direct calculations, we obtain that 
(6.5) 

G(<"'\<--i) ^ -(P(~^ sin2ax, ^ sin2^z), sin2nz, M^l^ sin27rz)*, 

2q; 2 2 2 

2 

G(^fi"'\Vn"'^) = -(P(^(cos2aa; + cos27rz),0),0,0)', 

2 

G(V'f^"'^^/'f^"^) = -(P(-^(cos2aa;-cos27rz),0),0,0)*, 

2a 

^/,/3iii,2 ,/3iii,2n /n^^^ • r, ■ n \ ^l(/3lll)7^ . „ ^2(/3lll)7r . j 

G(V'ii ,■011 ) = — (/'^(— sin2ax, — sin27rz), sin27rz, sin27rz) , 

2a 2 2 2 

and 
(6.6) 

G(<"-\ *fi--i) = -(P(-^ sin2ax, ^ sin2^.), ^ll^lii]! sin2^z, ^l^liill sin27rz)*, 

2a 2 2 2 

2 

G(V'fr''\ *fi"'^) = -(P(^(cos2ax + cos27rz),0), 0,0)*, 

2 

G(^/;f^"'^*f^"'^) = -(P(^(cos2ax-cos27rz),0),0,0)*, 

2a 

^/ ,/3iii,2 ,t,/3iii.2n /o/"^^ • r, ■ n \ C'l(/3lll)7r . G2(/3lll)7r . j 

G(V'ii , ^'ii ) = — (P( — sin2aa;, — sin27rz), sin27rz, sin27rz) . 

2a 2 2 2 

By dug) and (EH), we derive that for (j, Z) 7^ (0, 2), 

<G(<-''=\<-''=^),vI/f-^->ff=0, 

<G(^f/^<"-'=^),<i"''=^ >H=o, 
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where fci, fc2 — 1,2. Hence the first two terms in ()6.3p and (|6.4p are gone. Since the 
center manifold function contains only higher order terms 

*(a;iii, 2/111) = 0(|a;iiip, lymp), 

we derive that 



(6.8) 



< G(<i>, $), vl/^^i"^^ o(|xinr, 12/iiir), 



By dUll) and ([SH), only ym), (a^m, ym) and $^|(a;iii, ym) (where 

/?! = — = — 47r^ , (32 = = — 4T7r^, and f]^ — —<JJq2 = — 4cr7r^. ) contribute 

to the third order terms in evaluation of (|6.3p and (|6.4p . Direct calculations show 



< G«-\<--^),<^ >H= I' I' -^M§ii^sin^2.zdxdz = ^^^i^ii^, 

< G(<-'^ <-'^), <^ >H= Jj -^l^^sin^ 2.zd.dz = -^^^f^"^"' , 

< G•(<--^<--^),vI/(^^ >H= 0, < G(<"-\<"-2),vI/(^^ >ff= 0, 

< G«-\ <-'^), v,^^ >.= sin^ 2.zd.dz = ZME^^ 

< G(<--^<--^),vl/^^| 0, < G(<"-\<"-2),vi/^i >H= 0, 

< G(<"-\<"-i),vi/(^^ >H= 0, < G(<"-^<"-2),vI/(^^ >^= 0, 

< G(<--^<--^),vI/^^ >^= < G(<"-\<"-^),vi/g^ 



Applying Theorem 14.31 we obtain 



%2(a;iii,yiii) = — 777, (2^111 + Viii) + o{x^^^ +2/111), 

(6.9) ^/32^ ^ ^2(/3iii)7r 2 , 2 X , / 2 I 2 N 
^ ^ «'02(a;iii,yiii) = ^^^^ (a^iii +yiii) + o(a;iii +2/111), 

^02(^111,2/111) = o(a:;?ii +2/?ii)- 



By (I6?5])-(I6J]), we evaluate 
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11 'Mr'),^>H 
C27T' 



2a "2 
A2C2 



P2 



+yiii) + o(a;iii +^111) 



1 



8%/2 



/,/5iii,2 



*11 



11 



11 



>_ff = 



o(a;?ii +y?ii), 
+ Aid + A2C2), 



G(<"-\vI/?i--^),$(^^^o^2^ H-^^V-o^l >H +o(x?ii +y?ii) 



/3i 



-Gi^ 
2a 



^ AiGi 
'4a^ /5i 
1 



2a 

A2C2 
P2 



<'+o(a;^ii+ 2/^11 ) 

)(2^iii +yiii) + o(a;iii +2/111) 



8V2 



(AiGi + A2G2T-i)(x?ii + y2^i) + o(a;2ii + y^^J, 



< G($,^f{"'^),*^'- 



111 ^2 
11 



1 



o(a; 



2 

111 



yiii), 

2/lll): 



11 >H=^{-i + AiCi+A2C2), 



Plugging (|6.10p and (|6.1ip into (|6.3p and (|6.4[) respectively and applying Theo- 
rem 231 we get the reduced bifurcation equations: 



(6.12) 



dx 



111 



dt 



=/3iii(A)a;iii - i(5(A, 77)(a;?ii +Xiii2/iii) 
+ o(x?ii + y?ii) + 0(/3iii(A)(x?ii + yi^ii)), 



(6.13) 

where 
(6.14) 



dy 11 
dt 



=/?iii(A)2/iii - -5(A,7?)(x?ii?/iii +y?ii) 
+ o{xl^i + zj?ii) + 0(/3iii(A)(a;?ii + yJn)), 



5(A,r]) 



(A1G1+A2G2T-I) 



(3 + A1G1+A2G2) 
The following lemma determines the sign of 6{\^7f}. 
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Lemma 6.1. (1) Under the assumption ^57. 

f > if 71 < Tic, , 

"\ <0 tfri> T]c,. 
(2) // we replace r > 1 (t ^ a) in 15.1^1) . then 5{X,ri) > for all r/ > 0. 

Proof Step 1. Under the assumption A w = f + and /^m w 0. It 

suffices to prove the lemma for A = ^ + and /Jm = 0. Note that 
. crA J ^ ^ _ 

AiCi = -— ; ^"^^ ^2^2 — 7-3 ; 5^ ' 

(/3iii+7ii)^ +^7ii)^ 

Plugging A = ^ + and Pm — into the denominator and numerator of (5(A, rj) 
respectively yileds 

3 + AiCi + A2C2 « 3 + -i^^{a- - a^) 

>o, 

and AiCi + ^2^2X^1 « a7f/(l - r2)T-3(r;,^ - 



Hence the sign of S{X,r]) is determined by (ryci — 77). 
Step 2. If r > 1, we have 

3 + A,C\ + A2C2 « 3 + 7n (^^- - > 0, 

27 

and AiCi + A2C2T-1 » aj^^\—TT^ - r^{l - t^)t-^) > 0. 
Hence, S{X,ri) > for all 77 > 0. This completes the proof. □ 

7. Completion of the proofs 

We demonstrate the proof of Theorem 13.31 Theorem 13.51 can be proved in the 
same fashion. 

7.1. S'l-Attractor. First, Assertion (1) of Theorem \^ follows from (pT^ - lpT^ 
and Lemma [6Tl Then, by Theorem 14. 2i the equations bifurcate from (0, Ac) to an 
attractor Ea for A > Ac, which is homeomorphic to S^. 

7.2. Singularity Cycle. We shall prove that the bifurcated attractor 'Sx of (|2.4p - 
(|2.8p given in Theorem 13.31 is a cycle of steady state solutions. First, let 

H' = {([/, T, S)\u{^x, z) = -u{x, z)}, 

H[ = HiClH'. 

It is well-known that H' and H[ are invariant spaces for the operator + G given 
by (|2.10p in the sense that 

L\r^ + G : H[ ^ H' . 

It is clear that the first eigenvalue of L\^\H'i is simple when rj < rjc- By the 
Kransnoselski bifurcation theorem ( see aomng others Chow and Hale [2J and Niren- 
berg [lOp. when A crosses Ac, the equations bifurcate from the trivial solution to 
a steady state solution in H'. Therefore the attractor contains at least one 
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steady state solution. Secondly, it's easy to check that the equations ()2.4p - ()2.8p are 
translation invariant in x-direction. Hence if ?/;o(a;, z) — {U{x, z),T{x, z), S{x, z)) 
is a steady state solution, then ipo{x + p, z) are steady state solutions as well. By 
periodic conditions in x-direction, the set 

= {^o{x + p, z)\p e M} 

is a cycle homeophic to 5^ in Hi. Therefore the steady state of p.4p - (|2.8p generates 
a cycle of steady state solutions. Hence the bifurcated attractor T,\ contains at least 
a cycle of steady state solutions. 

7.3. Asymptotic structure of solutions. It's easy to see that for any initial 
value ijjQ = {Uq,To, Sq) S H, there is a time ti > such that the solution ip — 
{U{t,ipo),T{t,ipo), S{t,'i{jo j) is for t > ti, and is uniformly bounded in C-norm 
for any given r > 1. Hence, by Theorem 14.11 we have 

(7.1) lim min \\iP{t,iPo) - =0. 

We infer then from (|6.12p and (|6.14p that for any steady state solution cj) = 
(e, T, S) e Sa of (|2.4p - (|2.8p . the vector field e = (ei, 62) can be expressed as 



(7.2) 



{ei = -rV2sm{ax + 6) costtz + wi(a;iii, j/m, /Jm), 
62 = rcos(aa; + 6) sIuttz + W2(a;iii, j/m, /3iii), 

for some < 6* < 2tt. Here 



(7.3) 



^'i(a;iii,yiii,/3iii) = o(vAii(A)) for i = 1,2. 

Now we show that the vector field 
(7.4) eo = r\/2 sin ax COSTTZ, r cos aa; sin TTz^ 

is regular in = x (0, 1). The singular points of e are (x, z) = {{k + 5)\/2, 5), 
(/fcV2, 1) and (fc^/2,0) with fc e Z, and 

- \f2roL cos ax cos TTz y/2riT sin ax sin ttz^ 



detDeo{x, z) — det 



-ra sm ax sm TTZ ttt cos ax cos ttz 



r^TT^^O, for (x,z) = ((fc + i)V2,i), 

- r^TT^ ^ 0, for (x, z) = (fcV2, 0), (fc\/2, 1). 

Therefore the vector field (|7.4p is regular, and consequently, the vector field e in 
(|7.2p is regular for any Ac < A < Ac + e for some e small. It follows from Theorem 14. 71 
that the vector field e of (|7.2p is topologically equivalent to the vector field Cq given 
by (|7.4p . which has the topological structure as shown in Figure 



7.4. Proof of Theorem 3.5. Inferring from (|6.12p - (|6.14p and Lemma |6.1[ the 
equations bifurcate from (0, Ac) to a repeller for A < Ac, which is homeomorphic 
to S^. Since {U,T,S) = (0,0,0) is a global attractor for A near , there exists a 
saddle-node bifurcation point Ao in between and Ac- This completes the proof. 
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